NPS68-86-005 



NAVAL 



POSTGRADOAIE SCHOOL 

Monterey, California 




A NOTE ON PROCESSING SEA LEVEL DATA 
BY 

LAURENCE C. ^^REAKER 
JULY 1986 



Approved for public release; distribution unlimited 



Prepared for; Chief of Naval Research 
Arlington, VA 22217 



FedDocs 
D 208.14/2 
NPS-68-86- 



005 




NAVAL POSTGRADUATE SCHOOL 



iffoiTuU OfS-a-SO,-OOS 



Monterey, California 93943 



RAOM R.C. Austin 
Superintendent 



David A. Schrady 
Provost 



This report is for the research project "Ocean Prediction Through 
Observation, Modeling, and Analysis" sponsored by the Physical Oceanography 
Program of the Office of Naval Research under Program Element 61153N. 
Reproduction of all or part of this report is authorized. 



This report was prepared by: 



UNCLASSIFIED 






n i JD L iY MJOK LIDr^J - ^r 



REPORT DOCUMENTATION PAGE NAVAL POSTGRADUATE SCHOOi 

MONTEREY OA 30043-3 1 0 1 



la REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



'o RESTRICTIVE VIARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



2b DECLASSIFICATION /DOWNGRADING SCHEDULE 



3 DISTRIBUTION/ availability OF REPORT 

Approved for public release; 
distribution unlimited. 



4 PERFORMING ORGANIZATION REPORT NUMBER(S) 

NPS-68-86-005 



5 MONITORING ORGANIZATION REPORT NUMBER(S) 



6a NAME OF PERFORMING ORGANIZATION 

NAVPGSCOL 

Dept, of Oceanography 



6b OFFICE SYMBOL 
(If applicable) 

68 



7a NAME OF MONITORING ORGANIZATION 



6c. ADDRESS (City, State, and ZIP Code) 

Monterey, California 93943-5000 



7b ADDRESS (Gfy, Sfafe, and ZIP Code) 



8a. NAME OF FUNDING / SPONSORING 
ORGANIZATION 

Office of Naval Research 



8b OFFICE SYMBOL 
(If applicable) 

(1122 PO) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 

N000146WR24027 



8c. ADDRESS (C/fy, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 





PROGRAM 


PROJECT 


TASK 


WORK UNIT 




ELEMENT NO 


NO 


NO. 


ACCESSION NO 


Arlington, VA 22217 


61153N 


RR0310306 







n TITLE (Include Security Classification) 



A Note on Processing Sea Level Data 






12. £^ERS0NAL AUTHOR(S) 
T.anrp-nrp T . Rrpplrpr 








’3a. TYPE OF REPORT 
Technical 


13b TIME COVERED 
FROM Jul 86 TO Aug 86 


14 DATE OF REPORT {Yejr, Month. Day) 
1 August 1986 


15 PAGE COUNT 



’6 SUPPLEMENTARY NOTATION 



17 COSATI CODES 


F'ELD 


GROUP 


SUB-GROUP 















18 SUBJECT TERMS {Continue on reverse if necessary and identify by block number) 



'9 ABSTRACT (Continue on reverse if necessary and identify by block number) 

Procedures for processing hourly observations of tidal elevation from NOAA/NOS tide gauges 
for the study of subtidal variability are outlined. Processing procedures include an 
adjustment for datum reference, low-pass filtering to remove the diurnal and semidiurnal 
tides, and a static correction for atmospheric pressure. These procedures are illustrated 
using a sequence of tidal elevation data from Santa Monica, California for the period 15 
September to 15 December 1984. 



20 OlSTRlSUTlON/ AVAILABILITY OF ABSTRACT 

Eunclassified/unlimiteo □ same AS rpt □ otic users 


21 ABSTRACT SECURITY CLASSIFICATION 
UNCLASSIFIED 


22a NAME OF RESPONSIBLE INDIVIDUAL 
T..-mrpnrP. C ■ Breaker 


22b TELEPHONE (Inc/ude Area Code) 

(408) 646-3162 


22c. OFFICE symbol 
68Bv 



DD FORM 1473, 84 MAR 83 APR edition may be used until exhausted SECURITY CLASSIFICATION OF THIS PAGE 

All other editions are obsolete 



A G c 



:u Ri ' 



CuA^SlFlC* Tl'- 




A Note On Processing Sea Level Data 
L.C. Breaker 



Oceanography Department 
Naval Postgraduate School 
Monterey , CA 93943 

Abs tract . Procedures for processing hourly observations of tidal elevation 
from NOAA/NOS tide gauges for the study of subtidal variability are outlined. 
Processing procedures include an adjustment for datum reference, low-pass 
filtering to remove the diurnal and semidiurnal tides, and a static correction 
for atmospheric pressure. These procedures are illustrated using a sequence 
of tidal elevation data from Santa Monica, California for the period 15 
September to 15 December 1984. 
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Introduction 



The purpose of this note is to outline a set of procedures for processing 
hourly observations of tidal elevation to obtain information on subtidal 
vari abi 1 i ty . 



Recommended Processing Procedures 

Sea level data obtained from the National Ocean Survey of the National 
Oceanic and Atmospheric Administration (NOAA/NOS) tide gages is normally 
provided on an hourly basis with tidal elevations given in feet. The first 
correction to be applied adjusts the tidal elevations which are referenced to 
an arbitrary datum initially, to datums which are referenced to either Mean 
Higher High Water (MHHW), Mean High Water (MHW), Mean Sea Level (MSL), Mean 
Low Water (MLW), or Mean Lower Low Water (MLLW). NOAA/NOS provides the 
necessary data for each tide gauge location to make these corrections. The 
correction takes the following form; 

Adjusted Sea Level (feet) = Raw Observation (feet) - Correction (feet) 

To obtain MSL for example, 

MSL (feet) = Raw Observations (feet) - MSL Correction (feet) 

It is noted that because relative variations in sea level are generally of 
primary interest with regard to subtidal variability, no formal requirement to 
apply a datum correction in this case, exists. 

The second adjustment that must be applied removes the diurnal and 
semidiurnal tidal components. The Godin low-pass digital filter is used to 
remove these components (Godin, 1972). The response of this filter can be 
expressed as 
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where 



A • A /Z5-242 
25 24 



A = sin(rmAto) 
^ s 1 n(TiAaa ) 



i(TiAaay 

a = frequency in cycles per day 
At = time interval between observations in days 
Ap represents the arithmetric summation of n observations with the filter 
operating on the data three times, first using 25-point means (A^^ ) and then 
using 24-point means (A 2 ^), twice. According to Godin, the frequency band of 
this filter extends from 0 to 0.8 cycles per day (or 30 hours), beyond which 
the response is effectively zero. The response of this filter is shown in 
Figure 1. The half -ampl i tude cutoff is approximately 2.8 days (or ~66 hours). 
The amplitude of the diurnal and semidiurnal tidal components is reduced to 
less than Z% of the maximum value within the passband. A FORTRAN program 
implementing this tide-suppressing filter is contained in Appendix 1. 

A third adjustment must he made to sea level data to remove the effects 
of atmospheric pressure. According to Lisitzin (1974), in the stationary 
case, the sea surface is depressed by 1 cm for every increase in atmospheric 
pressure of 1.005 mbar (and vice versa). It is noted, however, that such a 
simple static adjustment for atmospheric pressure does not take into account 
the effects of surface wind. 

A relationship of the following form is recommended for making the sea 
level pressure (SLR) correction; 



SL(t)3^j = (SL(t)-^) - (SLP(t)-SLP) 



( 1 ) 



1.000 vice 1. OOOcm/1 . 005 mbar is used here 
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where 



SL(t)adj = Pressure corrected sea level (cm) 

SL(t) = filtered sea level (cm) 

^ = mean value of filtered sea level (cm) 

SLP(t) = filtered sea level pressure (mbar) 

SLP = mean value of filtered sea level pressure (mbar) 

The filtered sea level should first be converted to centimeters. Next, 
the SLP data should be low-pass filtered in a manner similar to that used for 
the sea level data. In this regard, it is recommended that the Godin filter 
be applied to the SLP data as well. Equivalent filtering of SLP has the 
advantage of producing a smooth time series with the same "time constant" or 
degree of smoothness as the sea level data. As indicated, both SL(t) and 
SLP(t) should be mean corrected in order to complete the correction for 
atmospheric pressure. Finally, whether or not the two terms on the R.H.S. of 
(1) are actually added or subtracted will depend on the sign convention used 
to define sea level and sea level pressure. The correct choice, of course, 
will produce an adjusted series with decreases (increases) in sea level that 
correspond to increases (decreases) in sea level pressure. 

An Example 

A 92-day sequence of hourly sea level data from the NOAA/NOS tide gauge 
at Santa Monica, California was acquired to illustrate the foregoing 
procedures. A datum correction furnished with the data from NOAA/NOS was 
applied initially to produce a sequence referenced to MSL. Next, the Godin 
low-pass filter was applied to the sea level data to remove the diurnal and 
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semidiurnal tidal components (Figure 2). The upper panel in Figure 2 shows 
the mixed character of this tidal record which clearly contains both diurnal 
and semidiurnal components. The lower panel shows the sea level record with 
the tides removed. To apply the correction for atmospheric pressure, hourly 
SLP data from a nearby MDBC environmental data buoy (NDBC Buoy 46011 off Pt. 
Conception) were acquired. These data were likewise smoothed using the Godin 
low-pass filter (Figure 3). The predominant diurnal component in SLP is 
effectively removed using this filter (lower trace). In the final figure 
(Figure 4), the mean-corrected SLP (upper panel) is subtracted from the 
mean-corrected sea level (middle panel) to yield sea level which has been 
corrected for atmospheric pressure (lower panel). 

Discussion 

Many different factors influence sea level. The time scales over which 
these factors are important vary from seconds to millenia. Only relatively 
short-term effects on sea level (< a few days) have been considered in the 
foregoing discussion. Although use of the Godin filter was recommended for 
removing the semidiurnal and diurnal tides, other tide-removing filters exist 
(Walters and Heston, 1982). In fact, Walters and Heston found that a 
so-called transform filter they describe generally had more desirable 
properties than the Godin filter. The only significant drawback they found in 
using the Godin filter is that in addition to removing tidal variations, it 
also attenuates variations in the 2 to 4 day range. For applications where 
information in this particular band is not of primary interest, however, the 
Godin filter is quite satisf actory . 

The simple adjustment for atmospheric pressure described here assumes 
that the direct response of the ocean to a change in SLP is hydrostatic. 
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Hamon (1966), however, has shown that such a simple one-to-minus one 
correspondence between SLP and sea level may not always be observed. Factors 

such as coastal trapped waves may contribute to the deviations from the 
hydrostatic approximation (Hamon, 1976). It is also possible to examine the 
phase relationship between SLP and sea level using the techniques of 
cross-spectral analysis (Hamon, 1966). When such an analysis is conducted, 
corrections for atmospheric pressure can be defined for a range of 
frequencies. Using this approach, Hamon (1966) found that changes in sea 
level often follow changes in SLP with lags ranging from 0 to 2 days. 
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FILTERED SEA LEVEL PRESSURE 
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Figure 4. Filtered sea level pressure (upper panel), filtered sea level 

(middle panel), and filtered sea level corrected for filtered sea 
level pressure (lower panel). 
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APPENDIX 1 



FILE: GCDIN 



F C R i y A N 



^ 1 



GODIN FILIEF APPLIES A 24 
CNCE TO PC'JPL'f Dii\. TPE 



TUAN THE FI-' 



POINT. 



HF ECXCAF I /.ICE AND A 25 H~ 
:?ST OHTPHT FCINI IS 3u rC: 



;f data fcints to 

POINTS 



R - 



DI!1ENSIO.N SL (3DID) 

INPUT ACTUAL IU:ibrP .. 

NPEAD =' N’J:13E5 C? INPUT DATA 
DO t)6 J=1.NFZAD 
EEAC(5,11) SL(J) 

6o CONTINUE 

C IN'S EFT ACTUAL DATA FO’^EAT IN 

11 fCa.1AT{DATA F CP. HA I) 

C FILTER THE DATA 

CALL BCXCAF(SL,NF£AD,2l,1,NCCT1) 

CALI 6C:<CA?(5L,NrUTl,24,1,NCUT2) 

CALL 3 CXC A F. ( S L , N C U T 2 , 2 5 , 1 . N C U T 3 
C WRITE OUT IllE FILIERcD DAiA 

WRIIE (0 . 14) 

14 FORH AT {1 X, ' NUHBSF OF FILTH "EC VALUES'//) 
WRITE(t),15) NCUT3 

15 FGEH A i (11 X ,15,//) 

CO o7 K=1,NCIJTJ 

4RITE(b,12) SL(X) 

67 CONTINUE 

12 FORHAI ( IX, F7. 1) 

STOP 

END 

SOdROUTINE boxcar (SL, NIN, LTH, IDEC, NOUT) 

CI.HENSICN SL(2500) 

NOUT = 0 
ISTAPT = 1 

20 lEND = ISTAPT LTH - 1 
IF (TEND. GI. NIN) GO TO 4 4 
3UH = 0. 

DO 3D I = I ST APT, IE ND 
SUE = S J H + S L ( I ) 

3D CONTINUE 

NCUT = 10 U I <• 1 

SL (NOUT) = S'M/FLCA: (ITil) 

ISTART = ISTAtPT + IDrC 
GO IC 20 
FEIUPN 
FNC 
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